Asevolution: a relativistic N-body implementation of the (a)symmetron

Abstract We present asevolution, a cosmological N -body code developed based on gevolution, which consistently solves for the (a)symmetron scalar field and metric potentials within weak-field approximation. In is dynamic can form non-linear structures. A cubic term added in symmetron potential to make symmetry-broken vacuum expectation values different, motivated by observational tensions late-time universe. To study effects of dynamics, we also implement constraint solver making use quasi-static approximation, provide options evaluating background evolution, including using full energy density averaged over simulation box Friedmann equation. The asevolution validated comparison with Newtonian ISIS that makes There found very small effect relativistic corrections our test simulations; it seen masses, not be accurately solved approximation; observe formation unstable domain walls demonstrate useful way identify them code. first consideration indicates are more asymmetron scenario.